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ABSTRACT 


Fe-Cr-Al  and  Fe-Cr-Mo  alloys  are  being  studied  for  their 
high  damping  characteristics  and  possible  utility  to  the 
Navy  in  the  noise  reduction  and  vibration  control  of 
shipboard  machinery  systems.  This  research  studied  the 
dependence  of  damping  in  these  alloys  on  heat  treatment  and 
cyclic  strain  amplitude.  Tensile  specimens  annealed  at 
900  C  and  1100  C  were  subjected  to  cyclic  axial  loading  at 
low  strain  levels  on  an  MTS  mechanical  testing  system, 
generating  load-displacement  hysteresis  loops  at  various 
frequencies.  These  hysteresis  loops  were  analyzed  to 
determine  the  damping  response  at  various  levels  of  loading. 
Damping  capacity  is  found  to  be  dependent  on  both  strain 
amplitude  and  frequency.  A  comparison  is  made  between  these 
results  and  previous  work  which  evaluated  the  strain 
dependence  of  damping  in  these  alloys  using  more 
conventional  experimental  methods. 
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I.  INTRODUCTION 


A.  GENERAL 

All  dynamic  mechanical  systems  produce  noise  and 
vibration  to  some  extent.  This  is  undesirable  in  a  number 
of  respects  Continuous,  large-amplitude  vibrations  can 
impair  the  performance  and  ultimately  shorten  the  service 
life  of  system  components.  High  noise  generation  from 
machinery  can  cause  hearing  loss  to  operating  personnel. 
Underwater  acoustic  emissions  from  the  propulsion  plant 
noise  of  Naval  vessels  may  increase  their  vulnerability  to 
detection  by  enemy  sensors.  For  all  of  these  reasons  the 
Navy  has  a  continuing  interest  in  all  methods  of  reducing 
the  noise  and  vibration  produced  by  its  shipboard  mechanical 
systems . 

Schetky  and  Perkins  [Ref.  1]  list  three  possible 
techniques  for  noise  and  vibration  control: 

1.  Isolate  mechanical  sources  from  surfaces  which  radiate 
transmitted  energy. 

2.  Dissipate  energy  within  the  mach:t  -ry  structure 
through  the  use  of  absorbing  pads,  ruboer  insulation, 
resilient  mounts,  etc. 

3 .  Attenuate  vibrational  energy  through  the  use  of 
energy-absorbing  (i.e.,  high  damping)  materials. 

The  first  two  methods  take  up  space  and  add  weight, 
respectively.  Thus  the  Navy  has  undertaken  research,  of 
which  the  present  research  is  a  part,  in  an  effort  to  design 
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“quiet"  alloys  which  possess,  along  with  the  required 
engineering  properties,  high  internal  damping 
characteristics . 

B.  BACKGROUND 

All  materials  will  damp  vibration  to  some  extent;  this 
is  a  property  reflecting  an  excited  component's  ability  to 
attenuate  vibrational  energy.  However,  most  common 
engineering  materials  have  very  low  damping,  capacities.  For 
example,  Reference  2  indicates  that  the  specific  damping 
capacity  (SDC)  of  most  steel  and  aluminum  alloys  suitable 
for  structural  applications  is  well  under  10  percent  at  room 
temperature.  Alloy  systems  in  which  high  damping  capacity 
(SDC  over  20  percent)  has  been  observed  include  those  in  the 
Cu-Mn,  Cu-Al-Ni,  Cu-Zn-Al,  Ti-Ni,  Fe-Cr-Al  and  Fe-Cr-Mo 
systems.  The  present  research  studied  certain  iron-chromium 
based  high  damping  alloys  in  which  the  Navy  expressed 
interest  in  terms  of  possible  application  to  shipboard 
mechanical  systems. 

Reference  3  presents  several  mathematical  models  for 
materials  undergoing  forced  vibration.  For  the  general  case 
of  a  sinusoidally  excited  structure,  the  Kelvin-Voigt 
single-degree  of  freedom  model  is  the  simplest  one  which  can 
still  be  considered  to  represent  a  real  mechanical 
situation.  The  differential  equation  for  this  model  is: 


mx  +  cx  +  kx  =  Fexp(iwt) 


(1) 


where : 


m 

c 

k 

x 

F 

Fexp(iwt) 

w 

•  •  • 

x  and  x 


mass 

viscous  damping  coefficient 
elastic  spring  constant,  or  stiffness 
displacement  of  the  system 
exciting  force  amplitude 

a  complex,  time-dependent  forcing  function 

excitation  frequency 

first  and  second  time  derivatives  of 
displacement. 


The  steady  state  solution  to  this  equation  is: 


x(t)  =  Xexp(i(wt  -  0)  =  Fexp(iwt)/ ( (k  -  mw2)  +  iwc)  (2) 


where: 


X  =  displacement  amplitude 

<j>  =  phase  angle  between  the  displacement  esponse 
and  the  input  excitation. 

The  viscous  damping  coefficient  can  be  represented  as: 

c  =  2mzw(n)  (3) 


where : 


z  =  damping  factor  ranging  between  zero  and  one 
(for  no  damping  z  =  0) 
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w(n)  =  the  natural  angular  frequency  of  oscillation 
(i.e,,  the  resonance  frequency  for  a  single 
degree  of  freedom  system) . 

This  simple  mathematical  model  is  considered  sufficient 
to  describe  most  of  the  physical  systems  described  in  this 
thesis. 

1.  Mathematical  Expressions  for  Damping 

Damping  capacity,  a  material’s  ability  to  absorb 
vibrational  energy  and  transform  it  into  other  forms  of 
energy,  can  be  measured  and  expressed  in  a  number  of  related 
quantities: 

a.  Logarithmic  decrement  (6)  for  exponential,  cyclic  free 
decay: 

6=  ln(Ai/Ai+1))  =  (l/n)  ln(A0/An)  =  2irz/  (1-z2)  V2  (4) 

where : 


n  =  number  of  cycles  between  A0  and  An. 

b.  Normalized  bandwidth:  Damping  measured  using  the  half 
power  point  method  for  evaluating  a  resonance  peak  in  the 
system's  frequency  response  to  excitation: 

Normalized  Bandwidth  =  (w2-w1)/w(n)  (5) 

where: 
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^v^vjvm 


a)l/W2  =  those  frequencies  framing  co  (n)  at  which 
the  stored  energy  is  half  its  maximum 
value  (see  Figure  (l)).  The  amplitude  at 

U1  an<^  u2  0*707  times  the  maximum 
value. 


c.  Quality  Factor:  the  quality  factor  Q  measures  the 
sharpness  of  a  resonance  peak  in  the  frequency  response: 


Q  =  w(n)/(ai2-u1)  =  l/2z 


d.  Internal  Fraction:  the  inverse  of  the  quality  factor, 
also  called  the  loss  factor: 


n  -  1/Q  -  (w2~W2)/u(n)  =  2z  (7) 

e.  Specific  Damping  Capacity:  (SDC)  represents  the 

energy  dissipated  per  cycle  for  an  oscillating  system  in 
free  decay: 


SDC  -  (Ai+12  -  A^/A^ 


(8) 


where: 

Ai+l'Ai  =  amplitudes  at  the  end  and  beginning 
of  ?  cycle. 

If  (Aj_+1  -  A^)  is  small,  the  SDC  can  be  approximated  as: 

SDC  =  2(Aj_+1  -  Ai)/Aj[  (9) 
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Normalized  Bandwidth  Used  in  the  Half-Power 
Point  Method  of  Measuring  Damping 
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SDC,  expressed  as  a  percentage,  is  related  to  the  other 
expressions  for  damping  as  follows: 

SDC(%)  =  200tr/Q  =  200  tt(  (w2-wi)/ w(n) ) 

=  2007rn  =  400TI6  (10) 

f.  Phase  Angle:  The  angle  by  which  a  system  response 

lags  behind  the  input  excitation: 

tan (4> )  =  (l/7r)ln(Ai/Ai+1)  =  <$/tt  =  1/Q  =  n  =  2z  (11) 

2.  Damping  Measurement  Techniques 

A  variety  of  techniques  have  been  developed  for 
measuring  the  damping  capacity  of  materials  and  structures. 
An  inverted  torsion  pendulum  apparatus  has  been  frequently 
used  [Refs.  4,5,6]  to  measure  the  log  decrement  of  a 
material  undergoing  free  vibration  (i.e.,  free  decay).  The 
single  cantilever  beam  resonance  dwell  apparatus  developed 
by  Bolt,  Beranek  and  Newman,  Inc.  [Ref.  7]  is  a  popular 
method  for  measuring  a  material's  damping  characteristics 
under  conditions  of  continuous  forced  vibration.  Previous 
research  at  NPS  on  high-damping  materials,  including  the  Pe¬ 
er  alloys,  have  used  a  modified  version  of  this  technique 
[Refs.  8,9,10,11].  A  brief  description  of  this  technique 
follows. 
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An  electromagnetic  shaker  drives  the  clamped  end  of 
a  cantilever  beam  specimen  with  a  continuous  forced 
excitation.  Accelerometers  mounted  at  the  root  and  tip  of 
the  beam  measure  the  response  of  the  beam  to  this 
excitation.  A  spectrum  analyzer  processes  the  signals  from 
the  input  (root)  and  output  (tip)  accelerometers,  producing 
the  frequency  response  of  the  beam.  The  half-power  point 
method  is  used  to  measure  the  damping  capacity  from  the 
resonance  peaks  which  appear  in  the  frequency  response.  A 
strain  gage  mounted  at  the  beam  root  measures  the  strain 
amplitude  experienced  by  the  beam  at  the  first  modal 
(resonant)  peak.  The  free  vibrational  decay  of  a  cantilever 
beam  can  also  be  used  to  obtain  the  log  decrement  of  a 
material  [Ref.  12]. 

The  generation  and  analysis  of  hysteresis  curves  to 
obtain  the  damping  capacity  of  a  material  will  be  described 
later. 

3 •  Fe-Cr-Based  Damping  Alloys 

Ternary  Fe-Cr-Mo  and  Fe-Cr-Al  high  damping  alloys 
have  been  developed  independently  by  Japanese  and  West 
German  companies.  This  study  considered  two  Fe-Cr-based 
alloys  similar  to  the  German  (Vacuumschmelze,  GMBH,  Hanau) 
Vacrosil-010  alloys.  Reference  should  be  made  to  Escue 
[Ref.  13]  and  Ahktar  [Ref.  14]  for  information  regarding  the 
corrosion  characteristics  of  these  alloys  in  seawater.  The 
chemical  composition  of  these  alloys  has  been  analyzed  by 
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two  independent  analysis  laboratories,  and  is  essentially  as 
follows  (numbers  listed  are  weight  percents) : 


TABLE  1 


CHEMICAL  COMPOSITION  OF  THE  ALLOYS  STUDIED 


'.  15] 


Element 


Fe-Cr-Mo 


Fe-Cr-Al 


Iron 


84.7 


85.0 


Chromium 


Molybdenum 

Aluminum 


Carbon 


11.7 


2.92 


0.007 


11.5 


2.95 


0.007 


Traces  (less  than  0.002%)  of  oxygen,  nitrogen,  sulphur, 
nickel  calcium  and  about  a  dozen  other  elements  were  also 
found  in  the  alloys.  Table  2  lists  some  representative 
engineering  properties  for  these  alloys: 


TABLE  2 

MECHANICAL  AND  ELECTRICAL  PROPERTIES  OF  THE 
ALLOYS  STUDIED  [REF.  16] 


Ultimate  Tensile  Strength 

Yield  Strength 

Density 

Young's  Modulus 
Electrical  Conductivity 
Relative  Permeability 


65,264  psi 
46,410  psi 

0.2422  lb/ in  (7.47  g/cm) 
24.38  x 106  psi 
0.96  x io6  mho/m 


.  i m  A.".  VAMif-lA 


The  Fe-Cr-Al  and  Fe-Cr-Mo  alloys  evaluated  here  are 
of  the  class  of  ferritic  stainless  steels.  Ferritic 
stainless  steels  have  been  less  widely  used  in  engineering 
applications  than  austenitic  and  martensitic  stainless 
steels  due  to  their  relatively  poor  weldability, 
susceptibility  to  embrittlement  and  notch  sensitivity. 
However,  they  do  possess  good  corrosion  properties,  and 
associated  with  their  ferromagnetism,  excellent  damping 
properties. 

At  room  temperature,  ferritic  stainless  steels  exist 
as  an  alpha  phase  solid  solution  with  a  body-centered  cubic 
(bcc)  structure.  Ferritic  stainless  steel  alloys  usually 
contain  very  little  dissolved  carbon  and  it  usually  appears 
in  the  microstructure  in  the  form  of  very  fine  chromium 
carbide  precipitates.  Carbon  makes  up  less  than  0.01%  of 
the  alloy  compositions  in  the  present  research.  For  a 
discussion  of  the  Fe-Cr  phase  diagram  as  it  pertains  to  high 
damping  ferritic  stainless  steel  alloys,  see  O'Toole  [Ref. 


The  high  damping  capacity  of  the  Fe-Cr  alloys 


studied  in  this  research  is  generally  accepted  to  be  a 
magnetomechanical  hysteresis  mechanism,  which  will  be 
discussed  later.  This  magnetomechanical  phenomenon,  and  the 
subsequent  damping  capacity  of  the  material,  is  influenced 


by  the  following  physical  parameters. 
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a.  Cold  Work  [Ref.  17] 


Damping  capacity  is  strongly  deteriorated  by 
cold  work.  A  reduction  of  more  than  five  percent  completely 
destroys  the  damping  effect?  however,  it  can  be  restored  by 
a  succeeding  heat  treatment. 

b.  Magnetic  Field  [Refs.  4,8,12,17] 

Damping  decreases  sharply  in  the  presence  of  a 
magnetic  field.  At  a  certain  field  lev^l  the  magnetic 
domains,  and  damping,  will  disappear  entirely.  Therefore, 
these  alloys  should  not  be  used  in  applications  where  they 
will  be  subjected  to  magnetic  fields  exceeding  5000  A/m. 
Experimental  studies  of  the  damping  response  of  these  alloys 
in  the  presence  of  varying  magnetic  fields  will  be  covered 
later. 

c.  Magnetic  (Curie)  Transformation  Point 

When  the  temperature  of  a  ferromagnetic  material 
is  increased,  the  added  energy  reduces  the  magnetic 
permeability  (magnetization)  and  permits  magnetic  domains  to 
become  randomly  oriented.  Above  the  magnetic  transition 
temperature  or  Curie  temperature,  the  material  becomes 
paramagnetic  (i.e.,  non-magnetic)  and  loses  all  magnetic- 
related  damping  properties.  The  Curie  temperature  is 
approximately  725  C  for  the  Fe-Cr-Mo  alloy,  and  is  above 
700  C  for  most  ferromagnetic  alloys. 
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d.  Heat  Treatment  [Refs.  5,8,9,17] 


Annealing  these  alloys  at  temperatures  in  the 
900  C  to  1200  c  range  improves  the  damping  capacity  over 
that  of  a  non-treated  alloy.  The  means  of  cooling  used 
following  solution  treatment  also  has  an  effect  on  the 
resulting  damping  capacity  of  the  material. 

e.  Stress  and  Strain  [Refs.  4,8,12,17] 

Numerous  experimental  studies  have  evaluated  the 
strain-dependence  of  damping  capacity  in  high  damping 
alloys.  What  follows  is  a  brief  summary  of  some  of  the  most 
significant  of  those  studies,  including  previous  work  at  the 
Naval  Postgraduate  School.  In  all  of  these  studies,  the 
general  behavior  these  materials  followed  was  an  initial 
damping  capacity  increase  with  increased  strain,  passing 
through  a  maximum  damping  at  a  certain  critical  strain 
amplitude,  and  then  decreasing. 

Kawabe,  et  al.  [Ref.  4]  investigated  the  strain 
dependence  of  damping  in  an  Fe-Cr-Al  alloy  ( S I LENT ALLOY , 
developed  by  the  Toshiba  Corporation) .  Using  the  inverted 
torsion  pendulum  technique,  damping  was  measured  over  a 
strain  range  of  0-300  microstrain.  Maximum  damping  was 
achieved  at  a  strain  amplitude  of  about  100  microstrain. 
Damping  measurements  were  done  in  the  presence  of  various 
magnetic  field  intensities:  0,  3020,  4540,  and  15,100 
Amperes  per  meter,  respectively.  Damping  decreased  with 
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increased  magnetic  field  intensity  for  all  strain 
amplitudes. 

Suzuki,  et  al.  [Ref.  12]  tested  the  damping 
capacity  of  sintered  Fe-Cr-Mo  alloys  by  measuring  the  log 
decrement  of  freely  decaying  oscillations  of  various 
cantilever  beam  specimens.  Experimental  results  were 
presented  in  figures  plotting  6  vs.  beam  amplitude,  alloy  % 
chromium,  sintering  temperature  and  relative  density. 
Additionally,  the  effect  of  an  external  magnetic  field  on 
log  decrement  values  was  measured  through  the  use  of 
solenoids  which  either  demagnetized  the  specimen  or 
magnetized  it  up  to  1591  A/m.  Damping  initially  decreased 
rapidly  with  increasing  magnetic  field,  reached  a  minimum  at 
about  4000  A/m  (at  50  to  60  percent  of  the  original  damping 
capacity) ,  then  increased  slightly,  levelling  out  at  the 
maximum  field  values  measured.  Damping  capacity  increased 
with  increased  deflection  amplitude,  passed  through  a 
maximum  value,  then  dropped  off.  A  strong  damping 
dependence  on  porosity,  percent  chromium  and  sintering 
temperature  were  also  noted  in  these  sintered  alloys. 

Schneider,  et  al.  [Ref.  17]  studied  the  damping 
properties  of  a  number  of  Fe-Cr  alloys.  SDC  vs.  stress 
amplitude  curves  were  obtained  using  the  log  decrement 
technique  for  the  Fe-Cr-Al  version  of  Vacrosil-010,  for  one 
hour  anneals  at  1000  and  1100  C.  A  maximum  SDC  of  40 
percent  was  obtained  with  the  1000  C  annealed  sample. 
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Damping  measurements  were  also  conducted  in  a  magnetic 
field;  it  was  shown  that  damping  in  Vacrosil-010  increases 
slightly  with  increased  field  intensity  up  to  2000  A/m,  then 
drops  off  sharply  once  the  field  is  large  enough  that  it 
fixes  the  magnetic  domain  walls  within  the  specimen. 

Masumoto,  et  al.  [Ref.  5]  did  an  extensive  study 
of  the  heat  treatment  and  strain  dependence  of  damping  in 
Fe-Cr-Mo  alloys.  A  number  of  different  alloys  were  studied, 
with  molybdenum  mass  percentage  ranging  from  four  to  10 
percent  and  chromium  mass  percentage  ranging  from  10  to  20 
percent.  One  hour  heat  treatments  were  done  at  900,  1000, 
1100  and  1200  c,  with  specimens  either  furnace  cooled  or 
water  quenched  following  heat  treatment.  Damping  capacity 
was  measured  using  the  inverted  torsion  pendulum  technique 
to  obtain  log  decrement.  Results  were  reported  by  plotting 
internal  friction  vs.  maximum  shear  strain  amplitude  over  a 
range  of  0-300  microstrain.  Damping  capacity  was  again  seen 
to  increase  with  increased  strain,  reach  a  maximum  and  then 
decrease  in  all  specimens.  For  the  alloy  most  closely 
resembling  the  mass  percentages  of  Vacrosil-010,  an  Fe-10  % 
Cr-4  %  alloy,  maximum  damping  was  obtained  at  about  80-120 
microstrain  for  all  heat  treatments.  For  furnace-cooled 
specimens,  damping  decreased  with  decreasing  annealing 
temperature,  with  a  large  drop  between  the  1100  C  and  1000  c 
curves;  maximum  damping  was  attained  using  the  1200  C 
solution  treatment.  Maximum  damping  in  the  water  quenched 
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specimens  was  obtained  at  1000  C,  followed  in  descending 
order  by  the  900  C,  1100  C,  and  1200  C  specimens.  Thus  a 
different  relationship  between  damping  and  annealing 
temperature  was  evident  for  the  two  different  cooling  rates 
after  the  anneal. 

The  next  few  paragraphs  summarize  the  results  of 
previous  research  at  the  Naval  Postgraduate  School  on  the 
heat  treatment  and  strain-amplitude  dependence  of  damping  in 
Fe-Cr-Mo  and  Fe-Cr-Al  high  damping  alloys. 

O'Toole  [Ref.  8]  tested  the  damping  capacity  of 
an  Fe-Cr-Mo  alloy  (similar  to  the  German  Vacrosil-010  alloy) 
using  the  modified  resonant  dwell  apparatus.  Cantilever 
beams  were  annealed  for  one  hour  at  800,  900,  1000,  and 
1000  C  and  either  furnace-cooled  or  water-quenched.  He 
obtained  damping  values  over  a  strain  range  of  0-250 
microstrain  for  all  specimens.  The  water-quenched  specimens 
displayed  SDC  values  in  the  5-15%  range  that  were  strain- 
independent.  Furnace-cooled  specimens  showed  a  strong 
strain-dependence  and  followed  the  previously  observed  trend 
of  increased  damping  with  increased  strain  until  a  maximum 
is  reached,  then  falling  off.  Table  3  summarizes,  the 
results  of  O'Toole's  damping  work.  O'Toole  speculated  that 
annealing  relieves  residual  stresses  which  would  reduce 
damping  in  non-annealed  specimens.  Furthermore,  the  quick 
cooling  rate  attained  through  water-quenching  reintroduces 
residual  stresses  in  the  material,  and  thus  explains  the  low 
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TABLE  3 


MAXIMUM  SDC  AND  PEAK  STRAIN  AMPLITUDE  FOR 
FE-CR-MO  ANNEALED  AT  VARIOUS  TEMPERATURES 
AND  FURNACE  COOLED  [REF.  9] 


Anneal  Temp.  (C) 
800 
900 
1000 
1100 


Max  SDC  (%) 
36.21 
49.01 
36.35 
54.72 


Peak  Strain  (*  10“6) 
155.0 
67.1 
126.4 
72.7 


damping  obtained  in  those  specimens.  The  drop  in  SDC  for 
the  1000  C  specimen  was  curious,  and  O'Toole  examined 
optical  micrographs  and  conducted  X-ray  diffraction  in 
determining  whether  a  phase  change  at  that  temperature 
contributed  to  the  decreased  damping.  No  evidence  of  a 
phase  change  was  found. 

Ferguson  [Ref.  9]  studied  the  relation  between 
heat  treatment,  strain  amplitude  and  damping  in  both  the  Fe- 
Cr-Mo  and  Fe-Cr-Al  alloys.  Using  a  different  resonant  dwell 
apparatus  than  O'Toole,  he  obtained  damping  data  over  a 
strain  range  of  0-50  microstrain  for  cantilever  beam 
specimens  annealed  at  950,  1000,  1050  and  1100  C.  Within 
this  strain  range  damping  values  increased  and  reached  a 
maximum  at  about  10  microstrain;  however,  only  for  the  Fe- 
Cr-Mo  alloy  annealed  at  1000  C  were  strain  amplitudes 
obtained  that  were  high  enough  to  display  the  post-maximum 
damping  decrease  found  in  the  other  studies.  He  obtained 


maximum  damping  values  in  the  SDC  range  of  45-60%  for  all  of 
his  specimens,  with  the  1050  C  annealed  specimens  yielding 
the  highest  damping  values  for  both  alloys. 

Ferguson  also  conducted  tensile  tests  on  these 
alloys  and  discovered  low  strain  nonlinearities  in  the 
stress-strain  diagrams.  He  attributed  these  nonlinearities 
to  the  hysteresis  phenomenon  inherent  in  these  materials. 
Ferguson  also  searched  for  evidence  of  a  phase  change  within 
the  range  of  the  annealing  temperatures,  using  dilatometry, 
with  negative  results. 

Cronauer  [Ref.  10]  used  the  same  resonant  dwell 
apparatus  as  Ferguson  in  testing  the  temperature  dependence 
of  damping  in  cantilever  beams  tested  at  temperature  ranges 
of  22  C,  36-41  C,  60-69  C,  and  92-109  C.  No  temperature 
dependence  of  damping  was  evident  over  this  range.  Cronauer 
concluded  that  the  damping  behavior  was  insensitive  to 
temperature  changes  up  to  109  C,  and  probably  up  to 
temperatures  approaching  the  Curie  temperature.  Maximum  SDC 
values  of  60-65%  were  obtained  at  a  peak  strain  amplitude  of 
five  microstrain. 

Patch  [Ref.  11]  developed  a  computerized  data 
acquisition  system  for  handling  damping  data  obtained  with 
the  modified  resonant  dwell  apparatus  used  by  Ferguson  and 
Cronauer.  He  also  attempted  to  image  the  ferromagnetic 
domains  via  transmission  electron  microscopy,  but  due  to  the 
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configuration  of  the  particular  TEM  used,  this  proved  t  be 
impossible. 


5 .  Microstructural  Damping  Mechanisms  in  Fe-Cr  Alloys 
According  to  de  Batist  [Ref.  18],  damping  in 
ferromagnetic  materials  occurs  primarily  as  a  result  of  the 
nonlinear,  magnetomechnical  hysteresis  phenomenon.  A 
secondary  possible  source  of  damping  is  the  presence  of 
dislocations  in  the  material.  Dislocations  in  a  crystalline 
material  may  be  expected  to  move  as  a  material  is  cyclically 
stressed,  resulting  in  internal  absorption  of  energy.  A 
third  possible  source  of  damping  is  due  to  the  effect  of 
phase  changes  in  the  material;  however  this  does  not  appear 
to  be  relevant  in  the  case  of  the  Fe-Cr-Mo  and  Fe-Cr-Al 
alloys  studied  here. 

A  brief,  straightforward  explanation  for  the 
basic  of  magnetomechanical  hysteresis  in  ferromagnetic 
materials  is  provided  by  Cochardt  [Ref.  19],  and  is  as 
follows; 

The  energy  dissipated  during  a  stress-strain  cycle  as 
a  result  of  the  magneto-mechanical  effect  is  generally 
caused  by  the  irreversible  magnetostrictive  strain.  Every 
ferromagnetic  material  consists  of  so-called  domains  which 
are  more  or  less  randomly  oriented  in  an  unmagnetized 
material.  However,  on  the  application  of  a  magnetic 
field,  or  a  stress,  these  domains  tend  to  align  themselves 
in  the  direction  of  the  field,  or  in  the  direction  of  the 
tension  strain.  This  movement  of  the  domains  results  in 
an  irreversible  change  of  the  dimensions  of  the  material 
which  is  called  'magnetostriction. '  Accordingly,  when  a 
stress  strain  curve  for  an  unmagnetized  ferromagnetic 
material  is  recorded,  more  strain  is  measured  than 
postulated  by  Hooke's  Law.  ...  On  relieving  the  load, 
the  elastic  strain  reverts  to  zero,  assuming  plastic 
strain  is  negligible,  but  the  magnetostrictive  strains 


18 


,  v^r  ■v  vit  «pvi  im  mjtJTi¥T.snirfi^.w\i 


«rflj^»raif;mJu\3oaju\jiJnmanc«Kiunarww*VJRnf«w 


remains  almost  constant;  in  other  words,  the  curve  follows 
a  hysteresis  loop.  The  greater  the  area  of  the  loop,  the 
larger  the  damping  capacity  of  the  material.  .  .  . 

Cochardt  expresses  damping  capacity  in  terms  of  the 
logarithmic  decrement  d  =  Aw/2w,  where  w  is  the  mean  elastic 
energy  of  the  specimen  and  Aw  is  the  energy  dissipated  in 
the  specimen  per  cycle  per  unit  volume.  He  also  presents 
the  concept  of  a  critical  stress,  at  which  damping  in  a 
ferromagnetic  material  reaches  a  maximum.  For  a  cantilever 
beam,  where  stress  conditions  closely  correspond  to  those  in 
vibrating  blades,  shafts,  etc.,  Cochardt’ s  development 
concludes  that  the  critical  stress  is  approximately  40 
percent  of  the  mF^.imum  normal  stress. 

Cochardt  also  observed  that,  due  to  the  magneto¬ 
mechanical  effect,  the  areas  of  magnetic  and  mechanical 
hysteresis  loops  are  proportional  to  each  other. 
Additionally,  he  noted  that  certain  high  damping  alloys  do 
not  exhibit  appreciable  damping  when  tested  in  a  strong 
magnetic  field.  Willertz  [Ref.  6]  observed  this  as  well. 
He  used  an  inverted  torsion  pendulum  to  measure  the  damping 
in  Fe-12%  cr  specimens  under  applied  axial  and  static 
torsional  stresses.  An  applied  magnetic  field  of  47,800 
Amperes  per  meter  was  found  to  completely  suppress  magneto¬ 
mechanical  damping  in  the  sample.  Subtracting  the  damping 
measured  in  the  presence  of  the  magnetic  field  from  that 
measured  without  the  magnetic  field  enabled  him  to  determine 
the  damping  solely  attributable  to  the  magneto-mechanical 
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phenomenon.  His  results  proved  that  magneto-mechanical 
damping  is  the  dominant  mode  of  damping  in  ferromagnetic 
alloys.  Work  by  Kawabe  [Ref.  4],  Suzuki  [Ref.  12]  and 
Schneider  [Ref.  17]  previously  mentioned  supported  these 
conclusions  and  provided  additional  experimental  work  on  the 
influence  of  an  external  magnetic  field  on  damping  in 
Vacrosil-010  and  similar  Fe-Cr-based  alloys. 

6.  Analysis  of  Stress-Strain  Hysteresis  Loops 

Nashif,  Jones  and  Henderson,  in  their  book  entitled 
vibration  Damping  [Ref.  20],  acknowledge  ferromagnetic 
hysteresis  as  being  the  dominant  mode  of  damping  for 
"typical  constructional  metal  alloys  and  even  some  high 
damping  alloys."  The  cusped,  nonlinear  shape  of  the 
hysteresis  curve  shown  in  Figure  2  is  cited  as  being  a 
common  feature  when  a  nonlinear  mechanism  is  responsible  for 
the  observed  damping.  In  the  case  of  so-called  linear 
viscous  damping  mechanisms,  the  observed  hysteresis  loop  is 
elliptical  with  a  straight-line  major  axis. 

Analytical  representations  of  a  nonlinear  hysteresis 
loop  have  been  developed,  many  of  which  are  presented  in 
Reference  20.  An  example  of  these  is  the  following: 

->* 

or  =  E{e  ±  v/n[e(o)-e)n  -  2n“1  e(o)]}  (12) 

or 
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Figure  2.  Typical  Nonlinear  Hysteresis  Loop  Obtained 
from  the  Stress-Strain  Response  of  a 
Ferro-Magnetic  Material 
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(13) 


a  =  E(e){e  ±  n(e)e(o)|l  -  (e/e(o))2|n) 

where : 

E(e)  =  the  strain  dependent  elastic  modulus 

a  =  stress  during  the  loading  part  of  the  cycle 
a  =  stress  during  the  unloading  part  of  the  cycle 
e  =  strain 

e(0)  =  maximum  strain 

n  =  an  empirical  hysteresis  parameter 
n(e)  =  the  strain  dependent  system  loss  factor 
v  =  Poisson's  ratio. 

In  plotting  a  vs.  e  for  all  instants  of  a  cycle  of  a 
sinusoidally  oscillating  system,  a  hysteresis  loop  is 
obtained.  The  energy  dissipation  during  a  cycle  of 
deformation  per  unit  volume  of  the  specimen  material  is: 

2tt/w  2 

Aw  =  <jx7de  =  /  (de/dt)  dt  =  irnEe(o)  (14) 

0 

The  maximum  energy  stored  in  the  specimen  is 

w  =  0 . 5Ee (o) 2  (15) 

Therefore  the  loss  factor  can  be  expressed  as 

n  =  Aw/27rw  (16) 
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and  the  specif  damping  capacity  is 


SDC  (%)  =  200  mi  =  100  (Aw/w)  (17) 

Nashif,  Jones,  and  Henderson  [Ref.  20]  give  many 
examples  of  how  hysteresis  loops  can  be  analyzed  to  obtain 
damping  parameters,  the  strain  dependent  elastic  modulus, 
etc.  The  loss  factor  and  modulus  are  dependent  on  both 
frequency  and  temperature;  however,  it  is  stated  that  and 
E  can  be  assumed  constant  over  a  limited  range  of 
frequencies  and  temperatures. 

Ray,  Wren  and  Kinra  [Ref.  21]  are  currently 
conducting  material  damping  measurements  on  an  Al  6061-Mg-Zr 
metal  matrix  composite  via  direct  quantitative  analysis  of 
hysteresis  loops  obtained  by  the  axial  loading  of  tensile 
specimens  in  an  MTS  mechanical  testing  system.  Their 
approach  utilizes  the  trapezoidal  area  approximation  in 
measuring  the  areas  on  a  stress-strain  diagram  within  and 
underneath  a  hysteresis  loop,  the  ratio  of  which  represents 
the  specific  damping  capacity  of  the  specimen  (see  Figure 
3)  .  Sinusoidal  oscillation  of  the  loads  produced  the 
observed  hysteresis  loops,  in  which  stress  and  strain  varied 
in  accordance  with  the  following: 

a  =  o(o)  +  o(a)  cos  wt  (18) 
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Figure  3.  Damping  Measurement  from  the  Areas  in  a 
Stress-Strain  Hysteresis  Curve 
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e  =  e(o)  +  e(a)  cos  (wt  -  $) 


(19) 


where: 

0(0), e(0)  =  the  mean  stress  or  strain 
0(a) ,e(a)  =  the  alternating  stress  or  strain 

0;e  =  the  instantaneous  stress  or  strain 

=  the  phase  lag  between  strain  and  stress,  and 
w  =  the  angular  frequency  of  oscillation. 

Their  specimens  were  tested  over  a  30-2000 
microstrain  strain  range  and  a  0.1-10.0  Hz  frequency  range. 
This  research  concentrated  on  determining  how  damping  varied 
with  the  matrix  orientation  within  the  specimen.  It 
demonstrated  the  usefulness  of  this  simple  but  reasonably 
accurate  method  for  measuring  material  damping  over  strain 
ranges  higher  than  those  previously  obtained  using  the  more 
common  log  decrement  or  half  power-point  methods  presented 
earlier. 

C.  OBJECTIVES 

Previous  damping  studies  at  the  Naval  Postgraduate 
School  have  used  the  resonant  dwell  apparatus  to  obtain  SDC 
as  a  function  of  cyclic  strain  amplitude  at  quite  low 
strains.  The  objective  of  the  present  research  was  to 
examine  other  techniques  for  measuring  the  damping  response 
of  high  damping  alloys,  to  use  higher  strain  amplitudes  than 
achieved  previously,  and  to  link  together  the  damping 
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results  obtained  for  all  the  ranges  of  strain  amplitude 
studied. 


II.  EXPERIMENTAL  PROCEDURE  AND  DESCRIPTION 

OF  DATA  ANALYSIS 


Button-end  tensile  specimens  were  prepared  in  jp  ceordance 
with  the  dimensions  and  specifications  indicct„-.d  in  Figure 
4.  Some  of  these  specimens  were  machined  from  an  ingot  of 
the  Fe-Cr-Mo  alloy;  the  others  were  machined  from  an  .ingot 
of  the  Fe-Cr-Al  alloy.  The  specimens  were  anne-xiwd  tc ohs 
hour  in  a  tube  furnace  under  flowing  argon  gas  (to  minimize 
oxidation  .f  the  specimen  surfaces  during  heat  treatment')  at 
temperatures  of  900  C  and  1100  C.  Specimens  were  allowed  to 
furnace  cool  to  room  temperature. 

An  MTS  (Mechanical  Test  System)  Machine  was  used  to 
subject  the  specimens  to  sinusoidally  alternating 
compressive  and  tensile  axial  loads.  Five  different  ioad 
levels  were  used,  corresponding  to  the  maximum  and  minimum 
stress  and  strain  levels  listed  in  Table  4.  The  highest 
stress  level  corresponds  to  about  38  percent  of  the 
approximate  yield  stress  for  these  alloys. 

Axial  displacement  of  the  specimens  in  response  to  the 
alternating  loads  was  measured  using  an  MTS  longitudinal 
extensometer  with  a  0.5  inch  gage  length.  Load-displacement 
data,  in  the  form  of  a  hysteresis  loop,  was  obtained 
graphically  using  a  plotting  device.  Additionally, 
numerical  load-displacement  data  were  obtained  using  the 
F_Collect  data  collection  program  developed  here  at  NPS  for 
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TABLE  4 


MAXIMUM  AND  MINIMUM 
USED  FOR  THE 

AXIAL  LOAD,  STRESS  AND  STRAIN  LEVELS 
GENERATION  OF  HYSTERESIS  LOOPS 

Max  Load  (lbs) 

Max  Stress  (psi) 

Max  Strain  (ps) 

400 

3636 

160 

800 

7273 

300 

1200 

10909 

525 

1600 

14545 

685 

2000 

18182 

850 

data  collection  on  the  MTS  testing  machine.  The  program  is 
run  on  a  Hewlett-Packard  9826  computer;  data  collection  is 
controlled  by  a  Hewlett-Packard  3497A  Data  Acquisition/ 
Control  Unit  and  a  Hewlett-Packard  3437A  system  voltmeter. 
Hysteresis  loops  were  generated  at  each  loading  level  for 
each  specimen  over  a  frequency  range  of  zero  to  ten  hertz. 
The  F_Collect  data  collection  program  is  listed  in  Appendix 
A. 

The  strain  energy  dissipated  per  cycle  (Aw)  and  total 
strain  energy  (w)  was  obtained  by  measuring  the  areas  within 
and  underneath  the  hysteresis  loops  using  the  trapezoidal 
area  approximation.  The  ratio  of  these  two  areas  ( av/w) 
represents  the  specific  damping  capacity  of  the  specimen 
under  those  loading  conditions  at  that  frequency.  Since  a 
ratio  is  calculated  the  load-displacement  hysteresis  loops 
can  be  used  directly  in  determining  damping,  as  all  units 
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cancel  and  the  SDC  ratio  is  dimensionless.  Maximum  strain 
values  per  cycle  were  obtained  using  the  maximum 
displacement  per  cycle  in  each  loop  and  the  strain 
relationship: 

e(true)  =  ln[(lQ  +  Al)/10]  =  Al/l0  =  e(eng)  (20) 

where: 

AL  =  instantaneous  displacement  at  maximum 
loading 

10  =  gage  length  =  0.5  inches. 

At  low  strain  levels  the  true  strain  and  engineering  strain 
are  approximately  equivalent.  The  maximum  strains  varied 
slightly  between  loops  obtained  at  different  frequencies  at 
a  given  loading;  the  maximum  strain  values  presented  in 
Table  4  here  are  averages  for  all  the  hysteresis  loops 
generated  at  that  loading. 
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A.  QUALITATIVE  ANALYSIS  OF  HYSTERESIS  LOOPS 

A  set  of  typical  hysteresis  loops  obtained  during  cyclic 
axial  loading  of  an  Fe-Cr-Al  tensile  specimen  over  a  10 
hertz  frequency  range  is  shown  in  Figures  5  and  6.  This 
particular  specimen  had  been  annealed  for  one  hour  at  1100  C 
and  furnace-cooled;  the  maximum  stress-strain  levels  were 
3620  psi  and  160  microstrain,  respectively.  This  figure  is 
representative  of  the  hysteresis  loops  generated  for  various 
specimens  at  various  loading  levels.  The  hysteresis  loop 
starts  out  very  narrow  at  low  frequencies,  representing  very 
little  energy  dissipated  per  cycle;  that  is,  low  damping. 
With  increased  frequency  the  loop  opens  up  and  becomes 
wider,  enclosing  a  greater  area  and  thus  indicating 
increased  damping.  At  a  certain  frequency  a  maximum  area  is 
obtained,  and  maximum  damping  is  achieved.  As  frequency  is 
increased  still  further,  the  loop  "collapses"  and  the  area 
enclosed  by  the  loop  is  reduced,  indicating  a  decrease  in 
damping. 

As  well  as  the  enclosed  area,  the  shape  of  the 
hysteresis  loop  changes  with  frequency.  At  the  lowest 
frequencies  examined  here  (1  to  4  Hz)  the  loop  has  the 
cusped,  nonlinear  shape  expected  with  a  material  in  which  a 
nonlinear  damping  mechanism  is  dominant.  At  slightly  higher 
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frequency  (5  to  10  Hz)  the  loops  become  elliptical  with  a 
straight-line  major  axis.  The  loops  in  Figures  5  and  6  also 
appear  to  be  getting  smaller,  but  this  is  actually  due  to 
the  plotting  device  being  unable  to  extend  through  the 
compressive  and  tensile  loading  maximums  at  high 
frequencies,  resulting  in  a  rounding  off  of  the  corners  of 
these  loops.  The  numerical  data  obtained  with  the  F_Collect 
data  acquisition  program  enabled  confirmation  that  the  high 
frequency  hysteresis  loops  were  elliptical  and  were  passing 
through  the  maximum  loads  and  displacements,  as  shown  for 
example  in  Figure  7,  which  is  a  representative  example  of  a 
high  frequency  hysteresis  loop  plotted  from  numerically 
acquired  data. 

As  previously  mentioned,  the  trapezoidal  area 
approximation  was  used  in  the  calculation  of  the  relevant 
areas  used  to  determine  damping.  This  method,  although 
tedious,  proved  to  be  a  fairly  accurate  method  provided  a 
reasonable  number  of  intervals  are  used.  Use  of  the  more 
accurate  Simpson’s  Rule  area  approximation  did  not  change 
the  calculated  damping  values  greatly. 

A  more  accurate  method  for  calculating  the  areas 
enclosed  by  the  hysteresis  loop  would  consist  of  a  computer 
program  which  would  use  the  numerical  data  obtained  by  the  F 
Collect  system  to  calculate  the  area  exactly.  such  a 
program  was  developed;  however,  problems  with  the  data 
acquisition  system  precluded  its  utilization.  The  main 
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DISPLACEMENT  (XI 0-4) 


Figure  7.  Load  vs.  Displacement  Hysteresis  Loop  from  a 
Fe-Cr-Mo  Tensile  Specimen  Annealed  for  One 
Hour  at  1100  C  Subjected  to  Cyclic  Axial 
Loading  of  ±  1600  lbs,  8.0  Hz 
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problem  with  the  data  acquisition  system  was  that  it 
typically  "dropped  out"  of  the  loop  after  one  quarter  or  one 
half  of  a  cycle,  and  generated  "garbage"  data  pairs 
thereafter.  The  data  acquisition  was  not  designed  to 
collect  load-displacement  data  pairs  at  the  speed  needed  to 
capture  a  hysteresis  loop.  As  a  result,  only  a  few  data 
runs  provided  a  complete  set  of  data  over  a  full  load- 
displacement  cycle.  Appendix  B  contains  the  area¬ 
calculating  program  used  in  conjunction  with  F_Collect. 

Since  the  numerical  data  were  incomplete  for  all  but  a 
few  instances,  and  since  the  high  frequency  hysteresis  loops 
were  not  accurately  plotted  by  the  plotting  device,  the 
areas  enclosed  by  these  loops  had  to  be  estimated.  This  was 
done  by  assuming  that  each  high-frequency  loop  was 
elliptical  in  shape  and  passed  through  the  load  and 
displacement  maximums  as  shown  in  Figure  7.  In  cases  where 
numerical  data  were  complete  and  the  actual  areas  could  be 
calculated,  the  estimated  areas  agreed  to  within  10  percent. 

Figures  8,  9,  10,  and  11  are  plots  of  specific  damping 
capacity  (SDC)  vs.  frequency  for  each  of  the  four  specimens 
tested.  For  each  specimen  the  damping  is  seen  to  be  greatly 
dependent  on  frequency,  increasing  rapidly  with  frequency 
until  a  maximum  is  reached  and  then  dropping  off  with 
increased  frequency.  The  frequency  at  which  the  maximum  SDC 
is  obtained  decreases  with  increased  loading;  i.e.,  the 
damping  vs.  frequency  curves  move  to  the  left  when  the 


Figure  8.  Specific  Damping  Capacity  vs.  Frequency  from  a  Fe-Cr-Al  Tensile 
Specimen  Annealed  for  One  Hour  at  900  C  Subjected  to  Cyclic 
Axial  Loadings  of  ±  400,  800,  1200,  1600  and  2000  lbs 
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Figure  11.  Specific  Damping  Capacity  vs.  Frequency  from  a  Fe-Cr- 
Specimen  Annealed  for  One  Hour  at  1100  C  Subjected  to 
Axial  Loadings  of  ±  400,  800,  1200,  1600  and  2000  lbs 


loading  level  is  increased.  Tabulated  data  on  which  Figures 
8  through  11  are  based  are  presented  in  Appendix  C. 

Maximum  SDC  was  obtained  at  the  1200  lb  load  level  for 
all  specimens.  This  corresponds  to  a  maximum  stress  and 
strain  level  of  10,909  psi  and  525  microstrain, 
respectively.  Figure  12  is  a  plot  of  maximum  SDC  vs. 
loading,  stress  and  strain  level.  This  plot  has  the  same 
shape  as  the  damping  vs.  strain  relationships  found  in  all 
the  previously  cited  experimental  results  measuring  damping 
in  high  damping  Fe-Cr-based  alloys.  Damping  is  seen  to 
increase  with  increased  strain  amplitude,  reach  a  maximum 
and  then  decrease. 

Table  5  lists  the  maximum  damping  capacities  obtained 
for  each  specimen  at  the  1200  lb  load  level. 


TABLE  5 

MAXIMUM  SDC  OBTAINED  BY  THE  FOUR  SPECIMENS 


Alloy 

Annealing  Temp. 

Maximum  SDC 

Fe-Cr-Al 

900  C 

64.36% 

Fe-Cr-Al 

1100  C 

62.24% 

Fe-Cr-Mo 

900  C 

54.22% 

Fe-Cr-Mo 

1100  C 

65.32% 

Due  to  the  small  number  of  samples  and  heat  treatments 
examined,  no  conclusions  regarding  heat  treatment  effects  on 
damping  can  be  made.  However,  the  dependence  of  damping  on 
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heat  treatment  has  been  firmly  established  by  previous 
research.  The  maximum  SDC  damping  values  obtained  are  all 
in  the  50-65%  range,  as  were  those  obtained  using  the  half 
power  point  method  for  cantilever  beams  randomly  vibrated  at 
much  lower  strain  amplitudes.  This  similar  damping  response 
within  completely  different  strain  ranges  (160-850  ys  here, 
compared  with  0-10  y  s  and  0-300  ys  in  previously  mentioned 
findings)  makes  it  difficult  to  link  the  independent  results 
in  order  to  make  any  general  observations. 

The  strong  frequency  dependence  of  the  damping  response 
of  these  materials  as  measured  using  the  hysteresis  loop 
analysis  was  unexpected  for  such  a  small  frequency  range. 
Concrete  physical  explanations  for  this  result  cannot  be 
readily  made,  although  it  is  suspected  the  frequency 
dependence  must  manifest  itself  on  a  sub-microstructural 
level,  influencing  the  interactions  of  atoms  and  the 
movement  of  the  magnetic  domains  within  the  material. 


Analysis  of  the  areas  defined  by  a  load  displacement 
hysteresis  loop  is  a  simple  and  reasonably  accurate  method 
for  measuring  the  damping  response  of  materials  at  higher 
strain  amplitudes  than  those  achievable  via  the  resonant 
dwell  or  inverted  torsion  pendulum  methods.  Further 
research  should  attempt  to  explain  the  frequency  dependence 
of  damping  measured  this  way  and  also  explore  possible 
reasons  why  the  damping  response  of  Fe-Cr  based  alloys  is 
similar  in  different  regions  of  strain. 

The  damping  peak  in  the  presented  SDC  vs.  frequency  data 
plots  is  very  sharp,  and  additional  hysteresis  measurements 
should  be  made  in  the  vicinity  of  the  maximum  to  define  the 
peak  damping  more  accurately.  Additional  hysteresis 
measurements  at  higher  and  lower  loading  levels  than  those 
done  here  could  be  done  to  extend  the  range  of  damping  vs. 
strain  data  obtained  using  this  method,  thus  defining  the 
SDC  vs.  strain  dependence  more  completely. 

Finally,  development  of  a  data  acquisition  system  able 
to  collect  reliable  numerical  hysteresis  data  from  the  MTS 
testing  machine  would  allow  exact  and  direct  calculation  of 
damping  capacity  from  the  hysteresis  loops. 


APPENDIX  A 


F  COLLECT  DATA  ACQUISITION  PROGRAM 


10  !This  program  is  F_COLLECT,  a  modified  version  of 
20  !a  program  written  by  Ken  Ellison  for  data 
30  ! collection  on  the  MTS  testing  machine 

40  Iwith  data  collection  controlled  by  the 
50  J HP  9826  computer,  HP  3497A  Data  acquisition/ 

60  ! control  unit,  and  the  HP  3437A  system 

70  ivoltmeter  (29  Sep  1987) 

80  ! 

90  !The  time  for  sampling  in  this  program  is 
loo  ! governed  by  the  total  number  of  samples 
110  ! based  on  2  bytes/point  and  the  delay 
120  !  The  sampling  is  uniform  so  by  sampling  2 
130  ! channels  the  time  between  data  POINTS  on 
140  isame  channel  is  2X  the  delay  time  and  the 
150  ! total  time  is  the  number  of  points  times 
160  !the  delay  time 

170  !  (500  points) x( . 002  sec  delay) =1.0  sec 
180  !  500  data  points=  250  data  pairs 

190  ! ******************************************** 

191  ON  ERROR  GOTO  Errcheck 

200  PRINT  "PROGRAM  NOW  LOADING" 

210  PRINT 

220  PRINT  "In  the  present  form  of  the  program  the" 

230  PRINT  "TOTAL  TEST  TIME  IS  ABOUT  1.0  SECONDS" 

240  PRINT 

250  COM  /Hp3054/  Sen, Dvm, Svm, Prt , Error , Err$ [ 6] , Ercnt$ [ 15 ] 

260  DIM  Y (1: 250) ,X(1:250) , In$ [1000] ,Out (1000) ,Name$ [50] 

270  BEEP 

280  INPUT  "What  is  your  load  range  (l-4):",Lr 

290  If  Lr=l  THEN  Lf=2  ! 2  kips/volt 

300  IF  Lr=3  THEN  Lf=l  !l  kip/volt 

310  If  Lr=3  THEN  Lf=.4  !.4  kps/volt 

320  IF  Lr=4  THEN  Lf=.2  !.3  kips/volt 

330  BEEP 

340  INFUT  "What  is  your  strain  range  (1-4):", Sr 

350  IF  SR=1  THEN  Sf=.01  ! Longitudinal  extensometer 

360  IF  Sr=2  THEN  Sf=.005 

370  IF  Sr=3  THEN  Sf=.002 

380  IF  Sr=4  THEN  Sf=.001 

390  PRINT  "Press  'continue'  to  initialize  the  STRAIN  &  LOAD" 
400  PRINT  "MTC  IN  HIGH  PRESSURE!!!" 

410  BEEP  300, .5 
420  BEEP  450, .75 
430  PAUSE 
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440  CLEAR  709  ! CLEARS  THE  DVM 

450  OUTPUT  724 ; "DCLAR"  i CLEARS  THE  DA/CU 

460  OUTPUT  724 ; "T1FN4SD. 2S" 

470  OUTPUT  709 ; "AC1AF4AL5AE2 " 

480  ENTER  724,*T,U,V 
490  J=0 

500  PRINTER  IS  706 
510 ! PRINT  T,U,V 
520  Disv=V 

530  DisO=Disv*SF  ! Initial  displacement 

540  Lodv=U 

550  LodO=Lodv*Lf  ! Initial  load 

560  PRINTER  IS  706 
570  PRINT 

580  PRINT  "Initial  displacement  is" ;DisO;"in  or" ;Disv; "volts 

590  PRINT  "Initial  load  is  ";Lod0;"kips  or  " ;Lodv; "volts 

600  PRINT 

610  PRINTER  IS  1 

620  BEEP  300, .5 

630  BEEP  450, . 75 

640  PRINT 

650  PRINT  "Press  ‘continue1  when  you  are  READY  TO  TEST." 

660  PRINT  "PRESS  'CONTINUE'  TO  START  DATA  ACQUISITION" 

670  PRINT  "THEN  PRESS  'HOLD*  ON  THE  MTS  TO  START  THE  TEST" 
680  PRINT 
730  CLEAR  709 

740  OUTPUT  709 ; "AC2AC1VNA" 

750  OUTPUT  724 ;»DCLF3D. OIOS" 

760  PAUSE 
770  ! 

780  !  Data  acquisition  portion 
790  1 

800  OUTPUT  7247 "T1N500SD. 002S" 

810  OUTPUT  709 ; "AF4AL5AE2 " 

820  ENTER  7247ln$ 

830  WAIT  1 
840  PRINT 

920  PRINT  "Data  acquisition  is  complete" 

930  PRINT  "Data  conversion  will  take  about  2  minutes 
940  ! input  from  the  DVM  the  array  of  data  strings 
950  J=o 
960  ! 

970  J  String  conversion  of  the  data  to  voltages 
980  ! 

990  •  *************ciiange  step  size  below  ****** 

1000  FOR  1=1  TO  1000  STEP  2  12  bytes/data  point 

1010  J=J+1 

1020  B1=NUM ( In$ [ I ] ) 

1030  B2=NUM ( In$ [ 1+1 ] ) 

1040  OUT  ( J)  =  (2*BIT  (Bl,  5)  -1)  *10A  (BIT(B1, 6)  *  (BIT  (Bl,  7)  -2)  ) 

1050  Out ( J) =Out ( J) *BIT ( Bl , 4 ) +BINAND (Bl , 15 ) + . 1*SHIFT (B2 , 4 ) 

+  .01*BINAND(B2,a5)  ) 
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1060  Out (J) =DROUND(Out (J) ,6) 

1070  NEXT  I 
1080  BEEP  330 , .50 
1090  ! 

1100  PRINT  '‘The  data  conversion  is  complete." 

1110  PRINT 
1120  1 

1130  INPUT  "Do  you  want  a  hard  or  soft  copy  (H  of  S)",Ans$ 
1140  PRINT 

1150  PRINT  "Is  there  any  interval  you  desire  to  see?" 

1160  PRINT  "Example:  every  5th  point  set" 

1170  PRINT  "enter  1  for  every  data  point,  5  for  every  5th" 
1180  PRINT 

1190  PRINT  "WARNING  -  THERE  ARE  260  DATA  SETS  TOTAL" 

1200  PRINT  "THAT  IS  A  LOT  OF  PAPER" 

1210  PRINT 

1220  PRINT  "PRINTOUT  IS  IN  GROUPOS  OF  50" 

1230  PRINT  "TO  PUNCHOUT  ENTER  '999'" 

1240  PRINT 

1250  BEEP 

1260  INPUT  0 

1270  If  Q<=0  THEN  Q=1 

1280  IF  Ans$="h"  OR  Ans$="H"  THEN 

1290  PRINTER  IS  706 

1300  ELSE 

1310  PRINTER  IS  1 
1320  END  IF 
1330  PRINT 

1340  PRINT 

1341  Punch=0 

1350  Fmtl:  IMAGE  DDDDD, 5X,MDD. DDDD, 5X,M. DDDDDDD 
1360  PRINT  "N", "LOAD", "DISP" 

1370  PRINT 

1380  FOR  1=1  TO  250 

1381  IF  Punch<>999  THEN 

1390  IF  I/50=INT (1/50)  THEN  INPUT  Punch 
1400  END  IF 

1410!****change  the  endng  value  above 

1420  X(I)=(Out(2*I) -Disv) *Sf  J Displacement 

1430  IF  I>1  THEN  X(I)=(X(I)+X(I-l)/2 

1440  Y (I) =Out (2*1-1) *Lf  ! Load 

1450  IF  I/Q=INT(I/Q)  OR  1=2  THEN 

1460  IF  Punch<>999  THEN  PRINT  USING  Fmtlyl, Y(I) ,X(I) 

1470  END  IF 
1480  NEXT  I 
1490  Dump;! 

1500  PRINTER  IS  1 
1510  ! 

1520  INPUT  "What  is  your  FILENAME? ",Name$ 

1530  PRINT 

1540  PRINT  "What  range  and  frequency  of  points  do  you  want 
to  save" 
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1550  PRINT  "'START,  STOP,  STEP"' 

1560  PRINT 

1570  INPUT  St,Sp,Stp 

1580  Siz=INT((Sp-St)/Stp)/15) 

1590  CREATE  BDAT  Name$,Siz 
1600  ASSIGN  @Path  TO  Name$ 

1610  FOR  J=St  TO  Sp  STEP  Stp 
1620  OUTPUT  @Path;Y(J)<X(J) 

1630  NEXT  J 
1640  PRINT 
1650  BEEP 
1660  PRINT 

1670  PRINT  "Tb'-.  program  is  now  complete" 
1680  BEEP 

1690  BEEP 

1691  GOTO  Compl 

1693  Errcheck: ! 

1694  If  ERRN=54  THEN 

1695  PRINT  "Duplicate  FILENAME-  TRY  AGAIN" 

1696  GOTO  1520 

1697  END  IF 

1698  Compl: ! 


APPENDIX  B 


AREA-CALCULATING  PROGRAM  FOR  USE  IN  CONJUNCTION 
WITH  THE  F  COLLECT  DATA  ACQUISITION  PROGRAM 


10  I********************************** 

20  !This  program  will  calculate  a  simple 
30  larea  under/between  the  curves  for 

40  ! hysteresis  loop  data  -  The  input  file 

50  ! should  be  in  the  form  LOAD, DISPLACEMENT 
60  !as  stored  by  the  data  collection  routine 
70  ! 'COLLECT'  or  ' F_C0LLECT ' 

80  i  11/5/87 

90  i************************************ 

100  !  Data  retrieval  portion 

110  INPUT  "What  is  your  INPUT  DAT  FILE  NAME" , Dat$ 
120  PRINT 

130  DIM  Lod(350) ,Displ(350) ,Strs(350) ,Strn(350) 
140  PRINT  "Is  your  data  disc  in?" 

150  PRINT  "Press  'CONTINUE*  when  ready" 

160  PAUSE 
170  L0=. 5 
180  Pie=3. 14159 
190  PRINT 

200  INPUT  "What  is  our  initial  diameter", DO 
210  PRINT 

220  AO=( Pie/4) *DOA2 

230  ASSIGN  @Pathl  TO  Dat$ 

240  FOR  N=1  TO  6000 

250  ON  END  @Pathl  GOTO  440 

260  ENTER  @Pathl ;Lod (N) , Displ (N) 

270  IF  Maxl<Lod(N)  THEN 
280  Maxl=Lod(N) 

290  Maxlpt=N 

300  END  IF 

310  IF  Minl>Lod(N)  THEN 
320  Minl=Lod(N) 

330  Minlpt=N 

340  END  IF 

350  IF  Maxd<Displ (N)  THEN 
360  Maxd=Displ (N) 

370  Maxdpt=N 

380  END  IF 

390  IF  Mind>Displ (N)  THEN 
400  Mind=Displ (N) 

410  Mindpt=N 

420  END  IF 
430  NEXT  N 


440  M=N-2 

450  ASSIGN  @Pathl  TO  * 

460  BEEP  300, 1 

470  PRINT  M;"  Data  points  have  been  read" 

480  PRINT  "Stresses  and  Strains  Are  now  being  calculated" 

481  PRINT  "PLEASE  WAIT" 

490  PRINT 

500  ! ***************************************** 

510  l CALC  OF  TRUE  STRESSES  AND  TRUE  STRAINS" 

520  ! 

530  FOR  N=1  TO  M 
540  Li=LO+Displ (N) 

550  Ai=(AO*L0)/Li  ! Based  on  const  vol  deformation 
560  IF  Li?LO  THEN 

570  Strn(N) =LOG(Li/L0)  !  In  tension 

580  ELSE 

590  Strn(N)=LOG(L0/Li)  J  In  compression 

600  END  IF 

610  Strs (NJ) =Lod (N) /Ai 
620  NEXT  N 
630  Area-0 
640  Totang=0 
650  Que=0 

660  ! *************************************** 

670  !  Area  under  the  stress/strain  curve  calc 
680  ! 

681  PRINT  "The  area  of  the  loop  is  now  being" 

682  PRINT  "calculated  -  PLEASE  WAIT" 

683  PRINT 

690  FOR  N=2  TO  M 

700  Angl=ATN(Lod(N-l)/Displ (N-l) ) 

710  Ang2=ATN (Lod(N) /Displ (N) ) 

730  IF  Totang<6. 283185  THEN 
740  Dstrn=Strn(N) -Strn(N-l) 

750  Avgstra= (Strs (N)+Strs (N-l) )/2 

760  Area=Area+ ( Ds trn* Avgs t r s j 

770  Que=N  !  The  last  point  for  1  full  loop 

780  END  IF 

790  NEXT  N 

800  !  Calc  of  the  area  fm  last  point  to  the  start  of  loop 
810  Dstrn=Strn(l) -Strn(Que) 

820  Avgstrs=(Strs(l)+Strs (Que) )/2 
830  Area=Area+(Dstrn*Avgstrs) 

840  ! ****************************** 

850  !  Printout  part 
860  ! 

870  PRINTER  IS  706 
880  PRINT 
890  PRINT 

900  PRINT  "For  the  hysteresis  loop  of  file  ";Dat$ 

910  PRINT  "The  total  area  enclosed  by  the  stress/strain  is:" 
920  PRINT  Area;"  ksi,in/in" 
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930  PRINT 
940  PRINT 

950  PRINT 

951  PRINTER  IS  1 
960  END 


APPENDIX  C 


TABULATED  SDC  AND  FREQUENCY  DATA 


1.  Specimen:  Fe-Cr-Al.  900 


Load:  ±  400  lbs 


FREQ  (HZ) 

SDC  (%) 

1-2.5 

8.20 

3.0 

11.87 

3.3 

36.09 

3.5 

42.20 

3.7 

46.19 

4.0 

47.92 

4.2 

29.01 

5.0 

27.54 

6.0 

19.71 

8.0 

13.31 

Load:  +  1200  lbs 


FREQ  (Hz) 

SDC  (%) 

1.0 

3.63 

1.5 

15.05 

1.8 

53.07 

2.0 

64.36 

2.5 

56.51 

3.0 

45.88 

4.0 

26.79 

6.0 

17.94 

8.0 

10.89 

Load: 

+  2000  lbs 

FREQ  (HZ) 

SDC  (%) 

0.5 

2.14 

1.0 

11.87 

1.2 

13.84 

1.5 

45.09 

1.8 

38.35 

(See  Figure  8) 


Load: 

±  aoo  lbs 

FREQ  (Hz) 

SDC  (%) 

1.0 

3.96 

1.5 

12.28 

1.8 

33.57 

2.0 

46.92 

2.3 

51.45 

2.5 

48.93 

2  s  9 

13.16 

4.0 

10.63 

6.0 

9.84 

8.0 

6.60 

Load:  ±  1600  lbs 


FREQ  (HZ) 

SDC  (%) 

1.0 

6.44 

1.3 

19.46 

1.5 

41.76 

1.8 

56.44 

2.0 

48.18 

2.5 

39.68 

3.0 

32.03 

5.0 

12.63 

7.0 

10.00 

FREQ  (HZ) 

SDC  (%) 

2.5 

39.61 

3.0 

21.94 

5.0 

11.94 

7.0 

8.04 
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Load:  ±400  lbs 


1100  C 


(See  Figure  9] 


FREQ  (Hz) 

SDC  (%) 

1.0-2. 5 

6.47 

3.0 

13.23 

3.2 

30.28 

3.5 

48.31 

3.8 

49.90 

4.0 

50.31 

5.0 

15.74 

8.0 

16.22 

Load: 

±  1200  lbs 

FREQ  (Hz) 

SDC  (%) 

1.0 

4.2 

1.5 

11.07 

1.8 

19.87 

2.0 

55.51 

2.2 

62.24 

2.5 

54.01 

3.0 

47.09 

4.0 

31.69 

6.0 

18.18 

8.0 

13.15 

Load: 

±2000  lbs 

FREQ  (Hz) 

SDC  (%) 

1.0 

6.36 

1.2 

29.03 

1.4 

35.22 

1.5 

40.94 

1.7 

35.65 

Load: 
FREQ  (Hz) 


1.0 

1.5 

1.8 
2.0 

2.5 

2.8 
3.0 
4.0 
8.0 


Load: 
FREQ  (Hz) 


1.0 

1.5 

1.6 

1.8 
2.0 

2.5 
3.0 
4.0 
6.0 
8.0 


FREQ  (HZ) 


2.0 

3.0 

4.0 

6.0 

8.0 


±  800  lbs 

SDC  (%) 

3.75 

12.05 

35.02 

51.46 

55.46 
51.32 
17.75 

9.94 

8.08 


±  1600  lbs 

SDC  (%) 

5.02 

14.71 

42.14 

52.58 
43.66 

42.58 
32.94 
32.94 
12.78 

7.90 


SDC  (%) 

30.97 

20.00 

13.42 

8.31 

5.38 


I 

1*1 

'»!« 


FREQ  (Hz) 


SDC  (%) 

3.20 

4.74 

32.48 

44.00 

44.98 

33.70 

28.83 

15.83 

10.70 


FREQ  (Hz) 


SDC  (%) 

2.84 

4.99 

47.88 

40.85 

31.75 

23.32 

18.95 

13.55 

12.46 


Load:  ±  2000  lbs 


FREQ  (Hz) 


SDC  (%) 


4.33 

5.42 

37.42 

35.68 

31.59 

31.62 

19.72 

11.07 


$ 


Specimen:  Fe-Cr-Mo.  1100  a  (see  Figure  li) 


Load:  ±  400  lbs 


ad:  ±  800  lbs 
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